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In this paper we study a curvaton model obtained by considering a probe anti-D3-brane with an- 
gular motion at the bottom of a KS throat with approximate isometries. We calculate the spectrum 
of curvature perturbations and the non-Gaussianities of this model. Specifically, we consider the 
limit of relativistic rotation of the curvaton brane which leads to a small sound speed, and thus it 
can be viewed as an implementation of the DBI-curvaton mechanism. We find that the primordial 
power spectrum is nearly scale-invariant while the non-Gaussianity of local type is sizable and that 
of equilateral type is usually large and negative. Moreover, we study both the theoretical and obser- 
vational constraints on this model, and find that there exists a sizable allowed region for the phase 
space of this model. 



I. INTRODUCTION 

The curvaton mechanism [l[ is an interesting proposal 
for explaining the observed scale-invariant primordial 
density perturbations in the framework of inflation (and 
some pioneering works on this mechanism have been con- 
sidered in Refs. @, I n this scenario, at least two 
fields are required, in which the curvaton field is subdom- 
inant during inflation but becomes dominant after the 
inflaton decays. When the curvaton is subdominant it 
provides entropy (iso-curvature) perturbations during in- 
flation, and afterwards, these entropy perturbations can 
be converted to curvature perturbations when the cur- 
vaton starts to dominate the universe. After the cur- 
vaton decays, the universe enters the standard thermal 
history, and then the primordial curvature perturbations 
lead to the formation of the large scale structure of our 
universe 0, @. 

Recently, a curvaton scenario realized in the frame of 
stringy inflationary models was presented in Ref. Q • As 
advocated in Ref. [|[ , a system of D- and anti-D-brancs 
in a warped throat can be made use of to provide an in- 
flationary stage in the early universe @, and its dynamics 
is described by a Dirac-Born-Infcld (DBI) action in the 
relativistic limit flol. [Tl| . Therefore, the curvaton mech- 
anism realized in this background is the so-called "DBI- 
curvaton" , where the curvaton field can be interpreted 
as a moduli parameter of a D-brane moving in a warped 
throat [l2l]. However, it is widely recognized that a brane 
inflation model of DBI-type often suffers from problems 
of backreaction both from the relativistic brane moving 
in the warped throa t FToL [T3 . [l4| and from the inflation- 
ary background (see [15| for an extended discussion), and 
so is often unable to provide a enough long inflationary 
stage required by cosmological observations [l|| [13, fl8j . 
The original DBI-curvaton model also suffers from the 
same problem as that of DBI inflation. 
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To circumvent this problem, one has to finely tune the 
precise shape of the potential of the inflaton in the pres- 
ence of supersymmetrically embedded D7-branes and an 
anti D3- brane localized at the tip of the warped conifold 
cone, and therefore the resulting model is rather delicate 
as studied in Refs. [TtJ H3, Hl| . Another approach to eas- 
ing this tension suggests that the motion of a probe brane 
along the angular coordinates in the internal space can be 
taken into account to prolong the inflationary period[22j 
which relaxed the theoretical constraints from the gravi- 
tational backreaction, and so extends the phase space of 
brane inflation [T3|. 

In this paper, we study a new DBI-type curvaton 
model established on the scenario of multiple brane in- 
flation proposed in Refs. [23l . |24| . We suggest that a 
probe anti-D3-brane with angular momentum on the tip 
of an approximately isometric warped throat can play 
the role of a curvaton field. Specifically, we consider the 
well-known Klebanov-Strassler (KS) throat [25| which is 
a nonsingular deformed conifold with its moduli being 
stabilized by fluxes [26l W\ . In this model, the anti-D3- 
brane always sits at the bottom of the deformed throat 
due to the sum of gravitational and gauge forces. Cor- 
respondingly its radial position is nearly fixed and the 
angular motion can be safely preserved during inflation. 

Since the specific realization of this curvaton mecha- 
nism involves at least two branes with the anti-D3-brane 
as the curvaton brane, this model can be viewed as an ap- 
plication of the multi- brane inflation model[H, [HI, |28| . 
and it naturally possesses some features of the multi- 
brane inflation model. For instance, the spectral index 
of iso-curvature perturbations generated by the curvaton 
brane during inflation has a slight red tilt deviating from 
scale-invariance, and is mainly dependent on the slow- 
roll parameter characterizing the variation of the Hub- 
ble parameter during inflation. At cubic order we find 
the nonlinear perturbations of a curvaton brane model 
have the combined features of a usual curvaton model 
and a DBI-inflation one. The non-Gaussianity of local 
type nicely agrees with the traditional result obtained in 
a general curvaton model; however, the non-Gaussianity 
of equilateral type is amplified by the small sound speed 
of the curvaton brane which is consistent with the usual 
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scenario of brane inflation. 

Our paper is organized as follows. In Section II we 
briefly review the dynamics of D- and anti-D-branes in 
the KS throat. In Section III, we present the realization 
of the DBI-curvaton in the KS throat and show that the 
curvaton brane can rotate in the relativistic limit. In 
Section IV we calculate the curvature perturbations and 
the non-Guassianities generated in this model. The con- 
sistency relationships are analyzed in Section V. Section 
VI contains discussion and conclusions. 



II. BRANE DYNAMICS IN A WARPED 
THROAT 

To begin with, we consider a flux compactification of 
type IIB string theory on an oricntifold of a Calabi-Yau 
three-fold . The geometry of this warped deformed 
conifold is viewed as a solution to type IIB supergravity 
with the presence of negative tension sources. It can be 
explicitly constructed by taking a stack of D3-branes with 
M units of the Ramond-Ramond (RR) fluxes F3 on the 
S3 (A-cycle) and K units of the Neveu-Schwarz-Neveu- 
Schwarz (NS-NS) fluxes H 3 on the dual cycle (B-cycle) 
which is the 5*2 times a circle extended along the radial 
direction. This gives, 
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F 3 = 2%M , 
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where a' is related to the string mass scale by a' = 1/ m 2 . 
When the fluxes are turned on, they warp the geometry 
of the conifold and the background metric of the 10D 
spacetime can be expressed as follows, 

ds 2 = tfg^dx^dx" + h- 2 (dl 2 + l 2 g mn dy m dy n ) , (2) 

where the warping factor h can be a function of the radius 
I and internal coordinates y m generically. Additionally, 
g^ u is the metric of 4-dimensional spacetime and g mn 
characterizes the internal space of the compact manifold. 

In the above system, we have only introduced the D3- 
branes. The geometry away from the tip of the conifold 
is an approximately Anti-de-Sitter (AdS) form with 
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and the characteristic length scale L is fixed by the flux 
number as 
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where g s is the string coupling which has to satisfy 
g s M ^> 1 to ensure the validity of the supergravity 
description 29] . In order to obtain a de-Sitter (dS) so- 
lution which may be used to realize inflation, we can in- 
troduce anti-D3-branes at the tip of the throat, which can 
lift the AdS vacuum such as the KKLT mechanism 12711. 



Note that, the warping factor at the infrared (IR) end of 
the throat takes 



hm — exp 
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and the IR cutoff of the radial coordinate lj 
{g s Ma')V 2 h IR . 

Now, we are interested in the dynamics of the branes 
in the KS throat. It has been shown that the D3-brane 
will sink to the tip of the throat, while the D3-brane is 
almost free except attracted by the Z?3-brane only. To 
have a quick look, we can write down the effective action 
for a D3-branc or an £>3-brane 130(1 . 



S = -T, 



det[G^} ± T 3 / d 4 a:(C4)oi23 , (6) 



with G^u being defined as 
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where + is for D3-brane, — for D3, and the X3 is the 

m 4 

tension which is given by T3 = m^hs g s - The equation of 
motion of the RR field yields, 

(C 4 )oi2 3 = V~9h* ■ (8) 

So, if the transverse fluctuations of the brane vanish, we 
will find that, for D3, these two contributions to the ac- 
tion cancel exactly; however, they could sink D3 to the 
bottom of the throat. 

As usual, if the KS throat is isometric, the warp factor 
h is independent of the angular coordinates y m . How- 
ever, in a general case, there are some corrections which 
make the warp factor dependent on angular coordinates. 
For example, since the compact C-Y manifold cannot 
have exact continuous isometries, the isometries of the 
bulk must be broken when the finite throat is glued on 
this bulk, and then the warp factor becomes angular- 
dependent. Since these effects are suppressed by powers 
of the warp factor, the correction to the warp factor can 
be described in terms of A(/i 4 ) ~ h^ R f(y) as analyzed in 
Refs. (U H2, [H| . However, we shall keep in mind that, 
A(/i 4 ) can not be too large, otherwise the geometry of 
the KS throat would be distorted seriously. 

A nonperturbative effect which stabilizes the Kahler 
moduli could bring a potential for the branes [HI, |35| |. 
The form of this potential depends on the precise embed- 
ding of the wrapped branes, and for embeddings which 
do not admit supersymmetric vacua on the tip, this po- 



tential could be angular-dependent. In Ref. [22[, the au- 
thors have provided an estimate of the nonperturbative 
effect which is dominated by the warp factor in the case 
of the Kuperstein embedding[36] of the D7-brane, where 
the angular degrees of freedom obtain a mass term, 
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in which fi measures the minimal radial location reached 

3/2 

by the D7-brane, and e ~ ^h IR characterizes the defor- 
mation of the coinfold. 

To combine the corrections mentioned above, we can 
obtain the effective action of the £>3-brane as follows, 

S=-T 3 J d 4 x^/-det[G^]-T 3 J d 4 xV=g(h 4 + V np ) , 

(10) 

where V np stands for the angular-dependent potential 
given by the nonperturbative effects, and so is a func- 
tion of the angular coordinates y m . Note that we have 
neglected the interactions among the curvaton and other 
fields for simplicity. In addition, even if the isometries of 
the throat are broken, the dynamics of the brane along 
the radial direction are almost the same as that in the 
isometric throat, unless the bulk effects can distort the 
throat geometry drastically. 

At the end of this section we would like to comment on 
the tadpole condition which is a limit on the net brane 
charge. Since the branes carry charges, they cannot be 
inserted into a compact space optionally. The extra brane 
charges must be compensated by the background, and 
thus the following condition has to be satisfied 

where \ is Ihe Euler number of the C-Y manifold, and 
Nd3 (-W753) stands for the number of the D3- brane (D3- 
brane). In the usual case, \ 1S °f order O(10 5 ) [37|, and 
so it puts an upper bound on the background charge 
number as MK < 10 5 . 

III. THE MODEL OF CURVATON BRANE 

Within the background introduced in the previous sec- 
tion, now we are able to construct the model of a DBI- 
curvaton from a Z?3-brane rotating in the KS throat. We 
consider a D3-brane with angular momentum at the tip 
of the KS throat which is approximately isometric. As 
shown in the previous section, this brane is localized in 
the radial direction. Note that in Ref. [14[ it was found 
that the angular momenta of the branes could be inflated 
away if these branes were used to drive inflation. In 
our construction, we assume that the inflationary back- 
ground is still driven by a D3-D7 system as introduced in 
Refs. [THHH (seejpUl] for earlier studies), the radial 
coordinate of the _D3-brane is nearly fixed, and there ex- 
ists an external potential from the nonperturbative effect 
of moduli stabilization which drives the curvaton brane 
to rotate. Consequently, provided the external potential 
persists sufficiently long, the angular momentum of the 
probe brane can safely survive during inflation. 1 In this 



case, the radial coordinate of the probe brane is around 
1 1 a, and thus we can study its dynamics along the angular 
directions at the tip of the KS throat and its application 
in the curvaton mechanism. The global scenario of our 
model is sketched in Fig. [TJ 




FIG. 1: A rotating D3-brane in a warped throat. In a de- 
formed KS throat, the branes are no longer pointlike in the 
extra dimensions. The circle stands for the angular motion of 
the brane at the infrared tip of the warped throat, Iir. 

Substituting the induced metric on the _D3-brane into 
its effective action (fT0|) . we can get a simplified form, 



S ~ T 3 I d 4 x^[-h 4 ^l + —(lj R g mn d^d^) 

-i^ + Vnp)}, (12) 

where we have neglected the higher derivative terms on 
spatial directions which are secondary on affecting the 
background evolution and the primordial perturbations 
at large scales. Note that, although the terms including 
higher spatial derivatives are negligible at linear and even 
cubic order, they could affect the nonlinear evolution of 
primordial perturbations at fourth order and so could 
leave a significant imprint on the trispectra of this model. 
This will be studied in following works in the near future. 

Moreover, there are three angular directions of the S3 
tip, and everyone is able to play the role of a curvaton 
field. We would like to pick up the flattest direction, and 
define it as the y coordinate, for simplicity. Consequently, 
we can define the curvaton field as, 

(j = liR\/%y. (13) 



This is different from Ref. [l4j ] in which their driving force is still 
almost radial and so leads to damping angular momenta. The 



non- vanishing angular degree of freedom in a warped throat was 
earlier applied to realize a bouncing universe 40, ill. 
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Finally, the effective action of the curvaton brane can be 
described by, 



-[h A (a)+V np (a)}\ . (14) 



As usual, we define a sound speed parameter of the cur- 
vaton field 



Cs = \ I 



■X 



(15) 



where X is define as X = —hip^crd^a). Further, we 
would like to define an effective potential of the curvaton 
field as follows, 

V{a) = T 3 (2/i 4 + V np ) , (16) 

and thus the corresponding effective Lagrangian can be 
written as, 



P = T 3 h 4 (l-c s )~V(a) , 



(17) 



of which the form is similar to that of the usual DBI 
model. 

We notice that a similar scenario involving a slow- 
rotating Z?3-brane was considered in Ref. [42| and its 
action was expanded to quadratic order under the as- 
sumption of the non-relativistic limit c s ~ 1. In the 
current work, however, we are going to consider the rel- 
ativistic dynamics of the model of the curvaton brane 
with the limit of c s ~ 0. In a flat Friedmann-Robertson- 
Walker (FRW) universe with the metric characterized by, 



d« 2 

as FRW 



-dt 2 + a 2 (t)S ij dx i dx i 



(18) 



the equation of motion for the curvaton field a is given 
by 



o + 2>H(t - —a-c s P a = 



(19) 



where the dot represents the cosmic time derivative, and 
P. a denotes the derivative of P with respect to the scalar 
field a. One may notice that, since the energy density 
of the curvaton is secondary during inflation, the hubble 
parameter H is determined by the background D3-D7 
inflationary model. 

Before studying the cosmological perturbations of this 
model, we would like to study the evolution of its back- 
ground equation under the relativistic limit. In the case 
of c s ~ 0, Eq. (HJ) yields 



(20) 



and it implies that the variation of the curvaton is 
roughly Act ~ —\{ 2 \fT 3 At for a fixed time interval At. 
Besides, as mentioned in Section 2, provided the geom- 
etry of the KS throat is not distorted too much, we 



still have the infrared warp factor hm. Following Refs. 
[23L |24| , we have an ansatz with the following form, 
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where we have set t — > —00 at the beginning of inflation. 
Moreover, Eq. (fT9"]l can be rewritten as 



+ v a = o 



(22) 

To apply the ansatz (|2"Tj) into Eq. (|2"2"|) , we then have 
the leading terms, which are the second term 



-3Hh^/% 



and the potential term 
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(24) 



respectively. The nonperturbative mass term is given 
by Eq. ([9]), as introduced in Section 2. The others are 
suppressed by the slow-varying factor h a and -hr (where 
\Ht\ ^> 1 is required by inflation). By matching the 
leading terms, we get the results p — 2 and a 
So we obtain the solution 
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(25) 



for a relativistic curvaton brane rotating in a warped 
throat. 

Notice that the value range of the curvaton is con- 
strained by the radius of the 3-circle at the tip. In the 
usual case, it restrict the variation of the coordinate along 
the angular direction Ay to be less than O(l). Conse- 
quently, to multiply the hubble parameter on the solu- 
tion (|25p . the angular restriction brings a constraint on 
the efolding number of inflation as follows, 



inf 



Hdt < H 



(26) 



In addition, the solution (f25|) gives a rolling-down be- 
havior for the curvaton field, which in the stringy frame 
shows that the rotation angle of the curvaton brane be- 
comes more and more parallel to the 3-cycle at the tip, 
and the angular-dependence will finally vanish. How- 
ever, the quantum fluctuation of the curvaton during in- 
flation may be in conflict with its classical variation so 
that hold the curvaton on the plateau of its potential 
longer. Therefore, with the help of quantum fluctua- 
tion there would be more angular modes preserved, but 
the model may suffer from the potential tension of large 
backreaction. 
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IV. PRIMORDIAL PERTURBATIONS OF THE 
CURVATON BRANE 



which implies that the power spectrum of the field fluc- 
tuations takes the form 



In the above section we have studied the background 
evolution of the curvaton brane during inflation. If our 
model indeed makes sense to the physics of the early 
universe, it is necessary to study the primordial pertur- 
bations generated by the curvaton field. This provides a 
potential discriminant between numerous curvaton mod- 
els in the literature. 



A. Quantum fluctuations during inflation 

To begin with, we split the curvaton field a into 
a(x,t) = a (t) + Sa(x,t), where a stands for the back- 
ground field which is homogeneous and isotropic, and 8a 
stands for the linear fluctuation which is caused by quan- 
tum effects. 

Following the formulae of Ref. [24J, we study the 
dynamics of quantum fluctuations in the spatially flat 
gauge, and assume no couplings to the background infla- 
ton. In this case, we have the equation of motion describ- 
ing the canonical perturbation variable in Fourier space 
as follows, 



8a 



(32) 



(c^-_)„ fc ~0 

z 



(27) 



where the prime denotes the derivative with respect to 
the conformal time r = J ^ . The canonical perturbation 
variable is defined as 



ee ac c 3/2 Sa k 



(28) 



and a background dependent function is introduced, 

z ee ac-^H^a , (29) 

in the above perturbation equation. 

One may notice that Eq (|2T[) has an asymptotic so- 
lution when we neglect the last term ^— , which implies 
\c s kr\ 3> 1, and it oscillates strongly like a sine function. 
This feature coincides with the slow roll condition de- 
fined above, which corresponds to the case in which the 
effective physical wavelength is deep inside the so-call 
sound horizon %. Therefore, the modes can be regarded 
as states in the adiabatic Minkowski vacuum during the 
sub-Hubble regime. Correspondingly, we can impose a 
suitable initial condition 



im 
v k 



-ic 8 kr 



(30) 



which corresponds to the Bunch-Davies vacuum. 

In the inflationary background, we have the following 
useful approximate relation, 



z 
z 



(31) 



after they exit the sound horizon. The subscript '*' 
stands for the sound horizon crossing time for the field 
fluctuation. 

Since during inflation the Hubble parameter is nearly 
constant and the field fluctuations are almost conserved 
after sound horizon crossing, we are able to calculate the 
spectral tilt of the primordial perturbations at the hori- 
zon crossing moment. This is given by 
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d\nP a 
G?ln k 



2d); — t— 3s* , 



(33) 



where we have defined a series of slow-variation parame- 
ters, 



H 

IP 



(34) 
(35) 



Note that we have ignored the square and the higher 
derivatives of the slow roll parameters in the above ex- 
pressions. The perturbation equation is derived under 
the assumption that the coupling between the adiabatic 
and isocurvature modes is negligible. A more detailed 
analysis of a similar scenario where an inflation model 
is based on multiple branes has been given in Ref. [24J. 
Thus we can naturally understand the difference in the 
results obtained in the current model and those obtained 
in [?J where multiple moduli degrees of freedom are in- 
volved in a single brane. Besides, to compare the r.h.s. 
of Eq. (|3"3")) with that in the usual curvaton model, we 
can see that the correction to the spectral index from the 
mass term disappears in our model, because it is strongly 
suppressed by a term proportional to . 



B. Curvature perturbations 

To proceed, we calculate the curvature perturbations 
generated by the curvaton field after inflation has ceased 
and the inflaton field has already decayed to radiation. 
During this period, the energy density of the universe is 
composed of p r which denotes the radiation energy den- 
sity and p a which is contributed by the curvaton. In this 
stage, the curvature perturbations are generated from 
the isocurvature modes since the pressure perturbations 
are non-adiabatic. This process ends when the perturba- 
tions become adiabatic again, which corresponds to the 
epoch of curvaton domination, or that of curvaton de- 
cay. The final curvature perturbations can be calculated 
at the moment H = T based on the assumption of per- 
turbatively instantaneous reheating of the curvaton. In 
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that case we can simply consider the component curva- 
ture perturbations Co- an d C7 on slices of uniform curva- 
ton density and radiation density separately. According 
to the assumptions of the curvaton mechanism, we can 
neglect the perturbations from radiation and thus have 
the curvature perturbation expressed as, 



c 



3(1 + w a p a ) 
4p 7 + 3(1 + w a )p a 



Ca 



(36) 



where w a = Pa/Pa is defined as the equation-of-state of 
the curvaton, and Ca is given by 



Ca = 



5 Pa 



3(1 + W a )p a 



(37) 



As emphasized in Section II, the nonperturbative ef- 
fect which stabilizes the Kahler moduli in the case of the 
Kuperstein embedding yields a quadratic potential for 
the curvaton brane as given in Eq. (|16|) . Therefore, after 
inflation the curvaton field arrives at its vacuum state 
and starts to oscillate when H = m np . At that moment, 
the sound speed c s goes to 1, and the energy density of 
the curvaton field redshifts proportional to a~ 3 . So we 
have an approximate equation-of-state w a = for the 
curvaton. In this case the energy density of the curva- 
ton is dominated by its potential, which in the oscillating 
epoch roughly takes the form 



Pa 



V(a) = 2T 3 h 4 IR 



1 



2 2 

m np a 



(38) 



where the first term is secondary in the oscillating epoch 
as will be explained in the next section. 

One can obtain the linear perturbation of the curvaton 
energy density in the oscillating epoch as follows, 



8 Pa 



ml p a Q 



z 8a a 



(39) 



where er osc is the amplitude at the beginning of its sinu- 
soidal oscillations. It can be viewed as a function of the 
curvaton field er» at Hubble exit, which is the only rele- 
vant quantity since the perturbations from the radiation 
fluid is supposed to be negligible. The curvature pertur- 
bation is calculated at the moment of curvaton decay, so 
it is convenient to define a parameter 



4p. 



7dec 



3p, 



(40) 



a dec 



which characterizes the fraction of curvaton component 
when it starts to decay. Consequently, at fixed p osc and 
Pdoc, we obtain the curvature perturbation at linear order 
as follows, 



37T(7osi 



(41) 



where we have applied the result presented in Eq. 
as well as the approximate relation cr osc CT> ~ 1 when the 
oscillation starts with H = m nv . 



C. Non-Gaussianities 

In the above section we have studied the primordial 
perturbations at linear order. It is necessary to ex- 
tend the theoretical framework of this model beyond 
the leading order. Particularly, the investigation of 
non-Gaussianities may provide a powerful discrimina- 
tor among numerous models describing the early uni- 
verse. For example, non-Gaussianities in a single canon- 
ical field inflation model were calculated in pioneering 
works [43|, |3] and found to be very small. Later, stud- 
ies of DBI inflation models (4f| predicted that a large 
non-Gaussianity of equilateral type can be realized in 
inflationary cosmology (46j. Recent developments in cos- 
mological perturbation theory shows that a sizable non- 
Gaussianity of local type can be obtained in bounce cos- 
mology but that it takes on negative values [47l. |48|. 



1. Local type 

To expand the perturbations to second order, we are 
able to investigate the features of non-Gaussianities in 
the curvaton brane scenario. At first glance, we consider 
the local form. The relative magnitude of the nonlinear 
perturbation is conventionally specified by a parameter 
/nl, which is defined by 



/i 



NL 



3 



c 2 



(42) 



In this case, it is convenient to use the 8M formalism 
to study the non-Gaussian perturbations by virtue of the 
separate universe assumption (4^ [5(| Hl|. Starting from 
an initial flat slice at time t osc , the curvature perturbation 
can be generally expanded as 



((t,x)=N jtT 5a+-Af tC 



(43) 



where Af(t, x) = ln[adec/flosc] is the amount of expansion 
to the final slice of uniform energy density at curvaton 
decay. To compare the first term in the r.h.s. ofEq. |43|) 
with the linear curvature perturbation obtained in Eq. 
(|4"Tj) . we can read N. a , — 2r a /3a osc . Differentiating this 
term with respect to cr* again and applying the definition 
°f /nLj w e can obtain 



NL 



6 Ml. 



5ry _5_ 
6 4r CT 



(44) 



This form precisely agrees with the result obtained in a 
normal curvaton modelfj], and more detailed studies of 
this parameter were present in [5(| H3, [53[ ■ 



2. Equilateral type 

The above result shows that the local non-Gaussianity 
in our model nicely coincides with that given in a general 
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curvaton model. For this case, we are unable to differ- 
entiate the model of curvaton brane from the standard 
curvaton model by observing the local non-Gaussianity. 
So we need to study more details of the nonlinear per- 
turbations generated in our model. In particular, we are 
interested in the non-Gaussianity of equilateral type. 

The presence of interactions in the perturbation La- 
grangian leads to non-Gaussianities. To start, we define 
the power spectrum and bispectrum as follows, 

<a,Cfc 2 ) = (2iT) 3 S 3 (k 1 + k 2 )P ( (k 1 ) , (45) 

(c fcl a 2 Gc3> = (2^) 3 < 5 3 (fc 1 + fc 2 + fc 3 ) 

xB c {h,k 2 ,k 3 ) , (46) 

and then these two spectra can be related in terms of a 
/c-dependent nonlinearity parameter Jnl, 



i>j/-vA,J,: ! ;-^(27r) 4 M^ /ArL (A :i ,A ;2 ,fc 3 ) , (47) 



in momentum space. 

It is easiest to work in the interaction picture, in which 
the three-point correlator to leading order is given by 



= iN% r J * dt'([5a kl 6a k2 6a k3 ,L 3 (t')}) , 

(48) 

where the SAf formalism was applied. For a generic cur- 
vaton model the integral stops at the moment of curva- 
ton decay. Here ti corresponds to the initial time be- 
fore which there are any non-Gaussianities. The square 
parentheses indicate the commutator, and L$ is the in- 
teraction Lagrangian which will be performed in the fol- 
lowing. In the second expression of Eq. (j4"8)) we have 
neglected the integral from the Hubble exit moment to 
the curvaton decay, due to the assumption of a curvaton 
model that the number of efolds from the end of inflation 
to the beginning of the oscillations is completely unper- 
turbed. 

If we have computed the three point correlators of cur- 
vaton field fluctuations, the nonlinearity parameter can 
be obtained by making use of the above equations. With 
an assumption of weak coupling between the fields, we 
perturb the sound speeds in the quadratic lagrangian and 
then obtain the lagrangian with the leading order terms 
up to cubic parts, 



2c, 5 cr' 



(49) 



Correspondingly, the dominant terms in the interaction 
Hamiltonian in Fourier space are given by 



Hi, 



dk 3 



a \(5& 3 + C -\k 2 5a5a 2 ) 



2c s 5 & 



(50) 



Then we decompose the field fluctuations in canonical 
quantization process, 



Scr k (t) = u{k)at + u* {—k)a/ - 



u(k) 



H 



:(1 + ic s kr)e 



— ic s k T 



(51) 



with the creation and annihilation operators defined by 
[a % ,at,] = d{k-k'). 

Since we have obtained the interaction Hamiltonian 
and the modes of the field fluctuations, we are now able 
to calculate the three point correlator of equilateral form, 



i l dt([dal, H int ] 



m 



-(2^) 3 <5(E^) 



(52) 



From Eqs. (fST))) and the above result, we can see that 
Hi n t ~ l/c s 2 . This roughly agrees with the result 
in usual single DBI inflation which is proportional to 
l/cs[ii]. To be explicit, we combine Eqs. ([¥T]). |47)l and 
(|48|) and eventually have the non-Gaussianity parameter 
of equilateral type, 



"H* 



NL 



2r, 



c 2 h 2 

' c s n IR 



(53) 



From the above analysis, we learn that the non- 
Gaussianity of equilateral type in the model of curvaton 
brane can be amplified both by small r a and c 2 s . We 
may understand this phenomenon as follows. The am- 
plification of non-Gaussianities from r a is inherited from 
the physical feature of a general curvaton model. More- 
over, when those perturbation modes exit the Hubble 
radius during inflation, the sound speed is also small due 
to the relativistic motion of the curvaton brane in the 
warped throat, and thus leads to significant equilateral 
non-Gaussianity signals. 

Note that in the above computations we assumed that 
the curvaton enters the oscillating epoch very soon after 
inflation. So an approximate relation o QSC (Jir ~ 1 can be 
used to obtain an estimate. Under this estimate there is 



5h 



lip 



(54) 



and so /] 



NL 



-5H*/2r a -c 2 mn P - It is worth extending the 



above analysis to more general cases. 



V. CONSISTENCY RELATIONSHIP AND 
CONSTRAINTS 

In this section, we shall discuss the consistency condi- 
tions and observational constraints on the phase space of 
the parameters appearing in our model. 

To our current knowledge, the most powerful observa- 
tional constraints on models of early universe are pro- 
vided by the CMB experiments. In order to link the 
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observations with the theoretical studies, let us make a 
short summary of the observable predictions obtained in 
our model. In our model, the amplitudes of the scalar 
and tensor perturbation spectra are given by 



p = r l m l P H * 



225-rr 2 hj R T 3 



Pt = 



2Hi 



(55) 



pi 



respectively, where m. p i — l/V87rG. The corresponding 
spectral indices are given by 



= 1 — 2e» 



3tS* 



n t = — 2e» 



(56) 



and the spectral tilt of the scalar perturbation spectrum 
is the same as that of the iso-curvature modes during 
inflation as shown in Eq. (|33[) . The tensor-to-scalar ratio 
takes the form, 



r T 



P T _ A50hj R T 3 



(57) 



l np'"'pl 



Going beyond the Gaussian statistics in the scalar per- 
turbation spectrum, the non-Gaussian fluctuations are 
characterized in terms of the nonlinearity parameter /nl- 
In the local and equilateral limits, the leading order con- 
tributions to this parameter are expressed as, 



/•local 
/NL 



5r g 
~6~ 



5 

4r„ 



pcquil 
J NL 



5if* 



2r a 



,(58) 



''up 



respectively. 

Since in our model the position of the curvaton brane 
is fixed at the tip of the warped throat, it suffers from 
very few theoretical constraints compared to usual brane 



54| 
56| 



inflation models, namely, the bounds from the throat 
and bulk volume [lU and the primordial tensor modes 
can be circumvented delicately. However, one still needs 
to be careful of the local gravitational backreaction of 
the probe brane on the throat. In the usual case this 
effect is negligible provided 4:ng s a' 2 -C Cgi 4 [l3|, and so 
gives a bound on the sound speed of the curvaton brane 
c s l > 0.6 /MK during inflation. 

Combining the Wilkinson Microwave Anisotropy 
Probe five year (WMAP5) data [57| with the distance 
measurements from the Type la supernovae (SN) [58[ 
and the Baryon Acoustic Oscillations (BAO) in the dis- 
tribution of galaxies the parameters of primordial 
cosmological perturbations have been accurately deter- 
mined by a group of stringent limits as follows: = 
(2.445 ± 0.096) x 10~ 9 , n c = 0.96 ± 0.013 at the pivot 
scale k = 0.002Mpc _1 , r T < 0.22 at the 2a level with a 
prior assumption of no running of the spectral index, and 
further -9 < < 111 and -151 < J^T' < 253 at the 

2<T level for the local and equilateral limits, respectively. 

Moreover, the number of efoldings of inflation is also 
constrained by cosmological observations. However, the 
reheating process after the primordial epoch is quite 
model-dependent and critically dependent on the reheat- 
ing temperature. Up to now, the allowed reheating tem- 
perature can be anywhere between lMeV and 10 16 GeV in 



principle. So it makes the measurement of the inflation- 
ary efolding number undetermined. In the current work, 
we take the simplest example of instantaneous curvaton 
reheating where the curvaton decays to standard-model 
particles momentarily once H = T. Therefore, we sim- 
ply assume the inflationary efolding number takes the 
value A/int = 50, which is a quite reasonable value al- 
lowed by observations. We would like to leave the more 
detailed study of the reheating process of the curvaton 
brane model to future studies. 

Additionally, there is one extra theoretical constraints 
from the requirement that the scalar perturbations gen- 
erated by the inflaton field should be negligible during 
the period of curvaton oscillations. We assume that the 
inflaton field satisfies the traditional slow-roll conditions, 
and hence the curvature perturbation during inflation are 

H 2 

given by the power spectrum P$ — » 2 . *i — ■ Similar to 

77 P t e * 

the definition of the tensor-to-scalar ratio, one can define 
a new parameter 



Ps 



rs 



225hj R T 3 



P C % r l m lp m li£* ' 



(59) 



which has to satisfy the condition rg < 1. 

Taking into account all the conditions introduced 
above, we find the available stringy parameters most 
stringently constrained by observations are g s and M. 
In Fig. O El and 2] , we present our numerical calcu- 
lations of the permitted regime for g s and M by fixing 
the total charge number of the background MK. Our 
results imply that the allowed phase space of the g s -M 
plane with a large MK is much bigger than that with a 
small one. This feature is similar to that of a traditional 
DBI inflation model, where a large value of MK can pro- 
vide a small enough warp factor; however, the limits on 
stringy parameters are much weaker in our case. This 
difference is mainly due to the particular feature of our 
model that the curvaton position keeps still during the 
epoch of inflation. 

In addition to the constraints on the stringy parame- 
ters performed above, we can also obtain the constraint 
on the slow-roll parameter as 0.0174 < < 0.0226. 
Moreover, in order to make sure that the equilateral non- 
Gaussianity generated in our model is consistent with the 
WMAP5 data, the fraction of curvaton density r a in the 
oscillating epoch needs to be of order 0(1). This is self- 
consistent with the previous assumption that the universe 
is dominated by the curvaton field after inflation. 

To conclude, in this section we have confronted our 
model with the latest cosmological observations. Under 
the priors of the inflationary efolding number and the 
background charge number, we derived that there is a 
sizable allowed region for the phase space of this model. 
It illustrates that cosmological observations can be nicely 
satisfied in our model. However, it is crucial to perform a 
global analysis which combines both the stringy parame- 
ters appearing in the model and the standard cosmolog- 
ical parameters without any artificial assumptions. We 



9 




10 



0.0 O.J 1.0 1.5 2.0 2.J 3.0 

Log ln M 

FIG. 2: Two dimensional constraints on g s and M from cur- 
rent observations, assuming a curvaton-dominated universe 
after inflation. The blue region is allowed by current observa- 
tions. In the numerical computation, MK — 10 3 and N — 50 
are chosen as priors. 




12 3 4 

Log ID M 



FIG. 4: Two dimensional constraints on g s and M from cur- 
rent observations, assuming a curvaton-dominated universe 
after inflation. The orange region is allowed by current obser- 
vations. In the numerical computation, MK — 2 x 10 4 and 
AT = 50 are chosen as priors. 



l.o I- h 




Log ln M 

FIG. 3: Two dimensional constraints on g a and M from cur- 
rent observations, assuming a curvaton-dominated universe 
after inflation. The green region is allowed by current ob- 
servations. In the numerical computation, MK — 10 4 and 
TV = 50 are chosen as priors. 



will leave this work to a future study. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper, we have presented a curvaton model ob- 
tained by considering a D3-brane in a KS throat with ap- 
proximate isometries from string theory. In this model, 
a D3-brane is stretched out along the 4D external space, 
and can be localized at the tip of the KS throat because 
of the gauge forces. Adding the nonperturbative effect of 
the Kahler moduli stabilization, its angular coordinates 
can play a role of the curvaton while the other degrees of 
freedom (for example, the radial motion of the Z?3-brane) 
are almost frozen. This scenario can be embedded into 
the framework of inflation models in terms of multiple 
branes, in which the Z?3-brane rotates relativistically. 

Based on the formulae developed in Ref . [23| and [HI , 
we have studied the background evolution and the pri- 
mordial perturbations of this curvaton model. The so- 
lution of the background dynamics of the curvaton is 
consistent with the assumptions of the nonperturbative 
potential of moduli stabilization, and so ensures the fea- 
sibility of our model. The process of converting the pri- 
mordial fluctuations of the curvaton into curvature per- 
turbations takes place after inflation, and is similar to 
what happens in the usual curvaton mechanism. We 
would like to emphasize that the final curvature per- 
turbations are calculated under the assumption that the 
curvaton decays immediately after the oscillation period 
ceases. Such a fast decay rate is also suggested by the 
stringy reheating configuration, in which one needs to 
turn on certain SUSY-breaking deformations of the KS 
background 60] in order to avoid redundant long-lived 
Kaluza-Klein (KK) modes[6l|,[63| and have a correct relic 
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density for DM candidates [60l [HI]. However, the pre- 
cise duration of the decay process is not addressed in 
the current work. This would alter the efficiency of the 
conversion of curvature perturbations and the precise de- 
termination of the inflationary cfolding number. 

Our model is related to many other interesting issues 
which deserve intensive studies in the future. For ex- 
ample, as mentioned below Eq. (|12p . this model is ex- 
pected to possess nontrivial kinetic couplings at fourth 
order from the geometric brane deformation. There- 
fore it may leave important signals on the trispectra 
which would render the model distinguishable from other 
multi-field inflation models with non-standard kinetic 
terms [si M, El E3, H M, M, [ill- In addition, a 
lesson from string theory suggests that a more reason- 
able expression of the multi-brane dynamics ought to be 
described b y a matrix in the framework of nonabelian 
background 72], which shows a potential connection be- 
tween the multi-brane scenario and the matrix inflation 



model investigated rece 

ntly[7l,[23- 
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